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We show that if a spectator linear isocurvature dark matter field degree of freedom has 
a constant mass through its entire evolution history, the maximum measurable isocurvature 
spectral index that is consistent with the current tensor-to-scalar ratio bound of about r < 
0.1 is about nj < 2.4, even if experiments can be sensitive to a 10^® contamination of the 
predominantly adiabatic power spectrum with an isocurvature power spectrum at the shortest 
observable length scales. Hence, any foreseeable future measurement of a blue isocurvature 
spectral index larger than ~ 2.4 may provide nontrivial evidence for dynamical degrees of 
freedom with time-dependent masses during inflation. The bound is not sensitive to the 
details of the reheating scenario and can be made mildly smaller if r is better constrained in 
the future. 


1. INTRODUCTION 

Although minimal single-field slow-roll inflationary seenarios [flUTOl ean sueoessfully provide 
a dynamieal explanation for the eurrently known features of the initial eonditions in elassieal eos- 
mologieal physies (e.g. the eosmie mierowave baekground (CMB) [|TTI - |25]| and large seale strue- 
ture it is natural to speeulate that more than one single real field is dynamieal during 

inflation. For sueh extra dynamieal degrees of freedom not to spoil the flatness of the inflaton 
potential, it is also natural to assume that they are very weakly eoupled to the inflaton (though this 
is obviously not a requirement). With this assumption, these extra dynamieal degrees of freedom 
behave as speetators as far as the inflationary dynamies is eoneemed. If one of these dynami¬ 
eal degrees of freedom is taken to be a weakly interaeting eold dark matter (CDM) field, then 
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there exists a well-known observable ealled the CDM-photon isoeurvature perturbations whieh 
beeomes observable (e.g. [|29] - l5n ) if the CDM field is suffieiently weakly interaeting and do not 
to thermalize. 

There are two broad eategories of sealar speetator field seenarios that ean produee observable 
CDM-photon isoeurvature perturbations: (i) linear speetators, sueh as axions [|52] - I54l . and (ii) 
gravitationally produeed superheavy dark matter seenarios, aka WIMPZILLAs [[55] - l60ll (for some 
reeent developments, see [lMT - [6^ f. Linear speetator fields are eharaeterized by having vaeuum 
expeetation values (VEVs) that are mueh larger than the amplitudes of their quantum fluetuations. 
The VEV oseillations generate the dark matter density in the universe today while the spatially 
inhomogeneous distribution of their energy-momentum tensors are determined by the quantum 
fluetuations. Sueh isoeurvature fluetuations are ealled linear beeause the the energy-momentum 
tensor inhomogeneity is approximately linear in the fluetuations, in eontrast with the the ease of 
gravitationally produeed superheavy dark matter seenarios. In this paper, we will foeus on the 
linear speetator seenarios and will drop the “linear” adjeetive.^ 

Seale-invariant isoeurvature perturbations with negligible eorrelations with eurvature perturba¬ 
tions are well eonstrained to be less than 3% of the adiabatie power [fT2l [TSl [T^ [65] - l70l . However, 
isoeurvature speetra with very blue speetral indiees ean be unobservably small on long wave¬ 
lengths, for whieh the measurements are strongly eonstraining, but have large amplitudes on short 
wavelengths, where the measurements are less eonstraining [1711 - 17^ . The ease of a blue speetrum 
is qualitatively different from a “bump” in the speetrum beeause bumps usually involve a red part 
as well as a blue part, and beeause the blue speetrum here is envisioned to have a qualitatively 
extended k-spaee range over whieh an approximately eonstant blue speetral index persists.^ 

One of the most natural models that ean produee large blue CDM-photon isoeurvature seenarios 
was given in iT/dll . This elass of models is eharaeterized by axions that have time-dependent 
masses due to the out-of-equilibrium nature of the Peeeei-Quinn (PQ) symmetry breaking field. 
Eor eonstant mass linear speetator fields, large blue-speetral indiees are diffieult to produee in 
observably large amplitudes beeause the energy density of the VEV dilutes away. An intuitive 
perspeetive is that the eloser the speetral index is to n/ = 1 (seale invariant), the more the field 
fluetuations behave like a frozen VEV, while the eloser the speetral index is to n/ = 4, the more 

* We briefly discuss what would happen with a quadratic isoeurvature scenario in the conclusions. 

^ Of course, from an observational point of view, this may not be easy to disentangle since observations have a finite 
k range. 
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the field fluetuations behave as partieles whieh ean be diluted away by inflation. 

Henee, a natural question, whieh is the subjeet of this paper, is what is the maximal measurable 
isoeurvature speetral index that ean be produeed by a eonstant mass speetator field in the eontext 
of slow-roll inflationary seenarios where the adiabatie perturbation speetrum originates from the 
inflaton field fluetuations. For a linear speetator seenario, we find that the maximum measurable 
speetral index in the foreseeable future is about n/ = 2.4 (where n/ = 1 eorresponds to seale in- 
varianee). Although measurability depends on the sensitivity of any given experiment, inflationary 
physios renders the dependenee of the experimental sensitivity to be logarithmie (to obtain some 
intuition, see e.g. Eq. ([50])). The bulk of this number originates from the ratio of the log of the dark 
matter density maximum enhanoement due to the dark matter diluting as a~^ (oompared to radia¬ 
tion diluting as a^‘^) and the number of efoldings neoessary for the inflationary seenario to explain 
the observed homogeneity and isotropy of the universe. A better eonstraint on the inflationary 
tensor perturbation amplitude r ean deerease this number, but the sensitivity is only logarithmie. 
If restrietions are plaeed on the maximum reheating temperature, then the maximum measurable 
speetral index also deereases. We will illustrate this by assuming a perturbative reheating seenario 
and assuming that the gravitationally suppressed nonrenormalizable operators of dimension 5 or 6 
are unavoidable. 

The number 2.4 is interesting beeause there are elaims in the literature [TT^ |75]-[73 that future 
experiments may be able to measure speetral indiees of nj >3. The results in this work demon¬ 
strate that if any of these experiments deteet a blue isoeurvature speetrum, then they may have 
uneovered evidenee for a dynamieal degree of freedom with a time-dependent mass. 

Before proeeeding, we note that the CDM-photon isoeurvature observable that we foeus on 
in this paper is distinet from the ^ eorrelator in the eontext of “heavy” masses diseussed e.g. in 
fTStiSOll and the ^-tensor eorrelators ll^ whieh in some eases ean also reeeive signatures from the 
isoeurvature degrees of freedom. On the other hand, these works all inelude the eommon theme of 
seeondary fields from inflation that ean leave a blue speetral eosmologieal observable signature. 

The order of presentation will be as follows. In See. we diseuss the eonstraints eonsidered 
in the speetral index maximization problem (there will turn out to be thirteen eonstraints). We 
then estimate the solution to the maximization problem analytieally in See.[^ Next, we solve the 
maximization problem numerieally in See. We then in See. [^give a brief review of why the 
axionie models that naturally have time dependent masses ean evade this bound and explain why 
this may be the most natural seenario to turn to if measurements are made of the speetral index 
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that are larger than ~ 2.4. Finally, we summarize and discuss caveats in the conclusions. 


2. MAXIMIZATION PROBLEM 

In this section, we define our class of models and the maximization problem at hand. In par¬ 
ticular, we provide a definition of a measurable blue isocurvature spectral index for a real scalar 
field X of constant mass that makes up a fraction of the total cold dark matter content through 
its background VEV oscillations, reminiscent of misaligned axion scenarios. 

We consider effectively single-field slow-roll inflationary scenarios, in which adiabatic cosmo¬ 
logical perturbations arise from the infiaton fluctuations. Here we define effectively single field 
to mean that a single field direction is important for the adiabatic inflationary observables. For 
example, hybrid inflation involves at least two fields, but during the slow-roll phase, only one field 
is dynamical as far as the adiabatic perturbations are concerned. 

In this context, consider a linear spectator isocurvature field X (see 11731 for a more precise 
definition) that is governed by the potential 

V(Z) = yZ", (1) 

in which m is a constant. Writing x = Zo(0 + the background equation of motion on the 

metric ds^ = dt^ — a^{t)\dx\^ is 

d^Xo + ^HdtXo + n?Xo = 0, (2) 

in which as usual H = d/ a. In accordance with the linear spectator definition, we assume that for 
the wave vector k in the range of isocurvature observable of interest, we have 

Xoitk) > (3) 

in which is the time when the mode k left the horizon (i.e. k = a{tk)H{tk)). The energy density 

in Xq oscillations that remains today is assumed to be part of the total cold dark matter content. We 
can then divide the dx (non-infiaton) perturbation into the adiabatic and non-adiabatic part in the 
Newtonian gauge as 5Xk = dx^^ + in which the nonadiabatic classical isocurvature field 

fluctuation 5^“^ = h^it) obeys the equation 

% -|- 3Hhi^ + m^h/^ = 0 


(4) 
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at the linearized level during inflation. If the approximate Bessel funetion solution index 
s/W — is real while the modes are subhorizon, then the square root of the ;t^-photon 

total isoeurvature amplitude is 


A?(t) 




Vejj^)2\h(t)\ 

lzo(<)l 


(5) 


(we have assumed the usual Buneh-Davies normalization \/^ in the limit oik/ (aH) —)■ 

oo), whieh remains frozen upon the horizon exit, where co^ is the eold dark matter fraetion eonsti- 
tuted by x, assuming all of dark matter is eold. More preeisely, the gauge invariant isoeurvature 
speetrum is (in the notation of [1731 1 


in whieh 


A2(k) = m;,2 


2'^-i 


2|r(v)| 


y/Tt ^ 

/ V 2:o(^^o) 

3 / 

4 m2 




f k\ 

UJ 


( 6 ) 

(7) 


in whieh Hitko) is the expansion rate when the ko mode leaves the horizon. Henee, for the blue 
speetral indiees that are of interest in this work, we have 


ni-\ 3-2v, 


( 8 ) 


in whieh v is a funetion of mass that also eontrols the time dependenee of the baekground field Xq- 
This elass of isoeurvature perturbations will be uneorrelated with the eurvature perturbations. 

We will eall nj — \ 0(1) > 0 a large blue speetral index, whieh eorresponds to m/H{t]^) ~ 

0(1). For the majority of this paper, we will take ko = h, whieh labels the longest wavelength 
mode relevant for CMB observations (around 0.002 Mpe^^), and we will assume SOe^t, -C 4 — «/. 
For brevity, we will also define Hi = 

As (0^ oc Xo{t) in Eq. ([^, and Xoi^) deeays exponentially during inflation whenever m/Hi ~ 
0(1), A^ ean easily beeome unmeasurably small for large blue speetral index seenarios. This 
suppression ean be partially offset by (k/ko)”^^^ enhaneements as long as 

eonstraint 1 : sj^{k^)/0)^ < 1 (9) 

to maintain perturbativity. In addition. Hi eannot in general be made arbitrarily large to make A^ 
large due to model building eonstraints sueh as the minimum number of efolds, reheating, and 
tensor perturbation limits. 

Given these eonstraints, a natural question arises: 
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• Given an experimental sensitivity parameterized by (whieh will be defined below), 
what is the maximum measurable «/ that ean be attributed to a eonstant mass speetator 
model in the context of effectively single-field inflation? 


This is the main question that will be answered in this paper, and the rest of the constraints (to¬ 
gether with Eq. Q) associated with maximizing u/ in Eq. Q will be laid out in this section. 

The main physics computation underlying this question is the determination of the time evo¬ 
lution of Zo(0 the time of reheating. The computation thus depends on the expansion rate 
H{t) during and after inflation. More specifically, the Xoi^) fi^ne dependence is governed by the 
time-coarse-grained amplitude of H because Eq. (|^ does not contain any derivatives of H{t). To 
cover a large class of slow-roll models economically (including both hybrid type and chaotic type), 
we consider a coarse-grained model space parameterized by Hi (the expansion rate when the the 
longest wavelength left the horizon), e^. (the potential slow-roll parameter when the longest mode 
left the horizon), and te — More precisely, we parameterize the expansion rate as 


H 


1 + I -HiSk. (4-4;)) 


( 10 ) 


which is continuous at tg? We will consider values that are consistent with the single-field 
adiabatic perturbation amplitude 


constraint 2: EkiiHi) = 


Ht 




( 11 ) 


in accordance with the spectator isocurvature paradigm considered in this paper. In the above, 
A^{ki) is the adiabatic spectral amplitude at the longest observable wavelengths, which we will 
take to be A|(k,) 2.4 x 10^^.^ As we seek a conservative upper bound on nj — 1, we will not 


^ Since we will never take the derivative of this function at 4 in the computation, the discontinuity of the derivative at 
f = 4 does not pose significant inaccuracies for the spectator field. This expansion rate fits the quadratic inflationary 
model to better than 10% during most of the time except at the inflationary exit transition where the fit degrades 
to 40% accuracy briefly at the transition point out of quasi-dS era. For inflationary models with smaller e*.., the fit 
is better, since this is a perturbative solution in Ekj- An alternative to this approach would be a numerical H time 
evolution sampling in the space of single held slow-roll models Il82ll87ll . We do not invest in the more numerically 
intensive approach since even an order 40% uncertainty in H amounts to an order 1% uncertainty in «/ — 1 in most 
of the parametric regime of interest. More discussion of this will be given later. 

^ This is consistent with current Planck measurements HU. A 10% change in this number only leads to less than a 
1% change in our results, while we are aiming for a 10% accuracy in «/ — 1. Hence, the precision of this number is 
not very important. 
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impose the adiabatic scalar spectral index constraint.^ The set of models that this parameterization 
excludes are those for which the quantities {£k., Hj, 4 } do not control Hitg), the expansion rate at 
the end of inflation. Such excluded models are somewhat atypical among known set of explicit 
effectively single-field models as they require new length scales (i.e. beyond Hi and tg) to enter 
the potential beyond those that are typically present in hybrid and chaotic inflationary scenarios. 
Furthermore, new length scales require yet another degree of fine tuning to fit smoothly with 
the t ~ tj^. time region where Eq. ( [TO] ) is guaranteed to be valid for effectively single-field slow- 
roll models. As we will discuss later, the maximum spectral index constraint does not sensitively 
depend on ei^.Hitg, which is fortunate since this parameterization is only 40% accurate for quadratic 
inflationary models near the time of the end of inflation. Note also that because we will impose 
the tensor perturbation phenomenological upper bound on Hj, the contribution to the spectral 
index will never be too big for phenomenological compatibility. 

In addition to the adiabatic constraint Eq. cn. we impose the inflationary condition that the 
number of efolds be larger than the minimum necessary for a successful cosmology: 


constraint 3: Ne= HiAlg 



>N 53 -h - In__ 

^3 lOiOGeV 


2 HejHute) 

3 IfliOGeV 


,( 12 ) 


in which 


Hg — Hi{\ Atg — tg 


(13) 


and we have taken the largest length scale to be kmin ~ IkHqqq. Note that in writing Eq. ( [12] ), 
we are neglecting contributions of order \n{chHi/{IHg)), in which c/, is an inflationary model 
dependent function of order unity. This leads to a systematic uncertainty with approximately a 2% 
error in the isocurvature spectral index bound. Note also Eq. ( [T^ is a non-linear constraint on Hi. 

We also impose the constraint that arises from assuming that there is at least one gravitational 
strength operator that can reheat the universe. Such assumptions are well motivated within string- 
motivated cosmologies (e.g. Il88] - l951l l and the weak gravity conjecture l!^ (for some recent devel¬ 
opments, see e.g. [l97] - l99l l. as well as generic expectations of interpreting gravity as an effective 
theory with the cutoff scale Mp. The minimum reheat temperature for a given Hg can be computed 
assuming a coherent oscillation perturbative reheating. Eor the inflaton field degree of freedom (p 


^ The imposition of the adiabatic spectral index constraint using a full chain of slow-roll parameter evolution sce¬ 
narios will not give a severe constraint on fg because of the large functional degree of freedom that exists in the 
inflationary slow-roll potential space, and its inclusion will obscure the presentation needlessly. 
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at the end of inflation to oseillate, we must have its mass nicp satisfy the eondition niq) > Hg. If the 
partiele deeay is through a dimension > 5 operator, then 






M] 


■2(«^-4) 


(14) 


is the gravitational deeay rate representing the “weakest” deeay rate where 5 is a phase spaee 
suppression faetor. For 2-body deeay, we expeets S ~ (8;r)^\ and we will take S as small as 
(0.1)^/(8;r) to get a eonservative bound. Sinee 


Trh 



(15) 


in whieh F is the total deeay rate, the bound F > and m,p > Hg lead to the following bound 


eonstraint 4: Hg < bound 



As we will see, for the maximal speetral index bounds at the highest reheating temperatures, this 
eonstraint is unimportant. A further eonstraint from reheating is that {Hi^Ek.jg} has to be ehosen 
for a fixed reheating temperature sueh that the energy at the end of inflation is large enough to give 
the total radiation energy: 

[3 

eonstraint 5: Trh < i — 1 y 

Here we have implieitly assumed Trh and Hg are sueh that eoherent oseillations of X oeeurs during 
the oseillation period of the inflaton. This eondition ean be written as 


3 

2 


HitRn) < m, 


(18) 


whieh ean be used to put a lower bound on the speetral index of 


eonstraint 6: 


n/ — 1 > 3 


\ 




(19) 


Although imposing eonstraint 6 seems artifieial sinee it is a simplifieation for ealeulational and pre¬ 
sentation purposes, the parametrie region where this bound is relevant is very similar to the para- 
metrie region where eonstraint 5 is relevant (i.e., it exeludes the similar {Hi^tg} region). Henee, 
there is no qualitative ehange in eomputing the maximum nj — Furthermore, we find in the 
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explicit numerical work that nj — I bound is lowered through constraint 6 by less than 1% which 
is below the systematic uncertainty in the computation. Hence, constraint 6 is a posteriori not 
important as long as constraint 5 is imposed. 

The absence of observed tensor perturbations yield the following phenomenological bound: 

constraint 7: Hi < Mp (20) 

in which is the bound on the tensor-to-scalar ratio (i.e. the ratio r < ~ IbCyt,). For the dark 

matter fraction to not exceed unity, we impose another phenomenological bound of 


constraint 8: (0^ < 1. 


( 21 ) 


We see that constraints 2 — 5 and 7 mainly arise from inflationary model-building consistency, 
while constraint 8 deals with dark matter phenomenology. 

We now turn to constraints on isocurvature perturbations in addition to constraints 1 and 6. 
Let us suppose that future experiments can detect isocurvature amplitudes above 

in which parameterizes the experimental sensitivity. Eq. (6) implies 


constraint 9: 


%2 




^|r(v)| \ r H,/i27t) \ 

J [ %o(fk,) 


ki 


i-v 


( 22 ) 


in which we have assumed 3/2 — v e^,.. We note that neglecting e^. in the spectral index is 
numerically valid to better than 2% level for the upper bound of interest. 

To see that constraint 9 controls the bound on the isocurvature nj — l that we are seeking, we 
note that if Xq oscillations occur before reheating, we have 

3 






(23) 


pR(Teq){Q.OM/{^b + ^DM)) ’ 
in which ELdm is the total dark matter fraction of the critical density today, is the total baryonic 
fraction today, and tgq is the time of matter-radiation equality. The prediction from the coherent 
oscillation perturbative reheating scenario takes the form 


co^=R 


2((^o)^)f=f«H 


m2 

IVAp 


m 


H{tRH) 


Trh 


leg 


in which 


^ _ ^b + ^DM 1 g*{TRH) g*s{Teq) 


3.94 1 

- ^ 0.23, 


^DM ^ g*{Teq) g^s{TRH) ^DM 3.38 6 


(24) 


( 25 ) 
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where g*{T) eounts the degrees of freedom in the radiation energy density pR, g*siT) eounts 
the degrees of freedom in the entropy density, and Tgq 0.8 eV is the matter-radiation equality 
temperature.^ As the solution to Eq. (|^ is approximately given by 


during inflation, while the next most important faetor is 


\ 2 
k 


-V 


"IZ 

= e 2 


'■"(0 


(27) 


with \n[k/ki] (number of efolds of inflation), we see that the magnitude of the left hand side of 
eonstraint 9 is eontrolled by CO^ and will be monotonieally deereasing as {nj — l)/2 inereases in 
the blue speetral parametrie region of interest.^ Henee, we eonelude that the maximum «/ — 1 is 
obtained when we saturate the inequality of eonstraint 9. 

It is also neeessary to eheek the eurrent phenomenologieal bound on the isoeurvature perturba¬ 
tions: 

- - (28) 


A2(ki) 




in whieh the eurrent phenomenologieal bound on for ki ^ 0.05 Mpe Ms ~ 0.2 at 95% eonfi- 
denee level ifT^ . Sinee A^(k) oc when eonstraint 9 is saturated Eq. (28) beeomes 




"izl 




< Eky. 


(29) 


To simplify the approximate phenomenologieal eonstraint parameterization, we ehoose k\ = kf. 

ki 


eonstraint 10: 


2 


< Ek:. 


Einally, we must also make sure we are in the linear speetator regime with our ehoiee of kn 

Elitk ) 

\ ^max / 


eonstraint 11: Zoih^J > 


2n 


(30) 


(31) 


and Xoihi) is not trans-Planekian: 


eonstraint 12: Xoihi) <Mp. 


(32) 


® Here we used gts{Teq) = 3.94 and g*{Teq) = 3.38. 

2 We see that intuitively when «/ — 1 = 0, the background held acts like a time independent constant while when 
nj—l —> 3^, the held behaves as a diluting gas of non-relativistic particles. 
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There is another uneertainty in constraint 9 that is associated with the fact that the Bessel func¬ 
tion mode functions are not obviously accurate solutions whenever the slow-roll parameter is not 
negligible. The limitations due to this issue were spelled out in iTTSl . A more accurate power-law 
expansion should have a fiducial value of = fcmax instead of ki. (The price that is paid for doing 
this is a complicated/numerical expression for Zo(^fcniax) terms of < Mp.) This will turn 

out to be an issue only for values of //, that saturate constraint 7 with rj, 3> 10”^ because 
evolves significantly in that case during inflation. To address this issue, for such worrisome situa¬ 
tions we therefore check the following constraint numerically 


constraint 9’: CO 






(33) 


^ \X0{tk,)\ ^-'^maxy-O 

involving a more accurate set of numerical solutions only. Finally, we note that constraint 9 also 
assumes that 


constraint 13: 


m 


H(tk 

\ f<-n 


^ 2 ’ 


(34) 


since only the non-decaying mode has been kept. We will see that in practice this does not pose a 
significant constraint. 

In summary, the problem of finding the maximally observable constant mass isocurvature spec¬ 
tral index nj for a given experimental sensitivity is to find the maximum nj that satisfies the 
constraints 1-13 given above. 


3. ANALYTIC ESTIMATE 


In this section, we provide an analytical estimate of the solution to the n/ — 1 extremization 
problem presented in Sec. We begin in Section [3T] by giving a crude estimate of the maxi¬ 
mization problem that is obtained by neglecting the slow-roll parameter e^t,- In Section 13.21 we 


then obtain an analytic perspective of the effect of turning on the slow-roll evolution of H and 
the non-linearities of the problem. For example, we will see that the Hi may not quite saturate 


constraint 7 for the largest spectral index, in contrast with the estimate given in Section 3.1 and 
this turns out to be significant for the accuracy of the approximation of the spectral index used in 
constraint 9. Sec. will involve a numerical solution to the constrained maximization problem 
without resorting to the analytic arguments presented in this section. 
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3.1. Without Slow-roll Evolution 


As previously discussed, the maximal spectral index results when constraint 9 is saturated. To 
evaluate constraint 9, we need to determine {{Xo) )t=tRH- this section, we will estimate this 
quantity to obtain a qualitative understanding of the parameters involved. 

Let us neglect the slow-roll evolution of H and assume that Xo coherently oscillates just at the 
end of inflation. We can then estimate 


{{Xof)t=tRH-C 


_ Hite) _ 


c = 


-{ni-\)Ne], 


(35) 

(36) 


in which Ne ~ Hi(te — t/t,) is the number of efolds of inflation. Through standard cosmological 
scaling, this yields the dark matter fraction to be 


co^ ~ R- 


xiitki)oxp[-{ni-l)Ne] 


Ml 


m 

2 

_ 1 


1 

\TgJ 


(37) 


Now, noting that m///,- ~ 0(1), and that the greatest n/ — 1 sensitivity comes from the exponential, 
we find (assuming constraint 1 is satisfied) 

( 




55 


2N, 


ki 


l + lln 


in-s? / QrM ) 

J f Trh \ 




\ lOiOGeV ) 




_ '^mm 


(38) 


in which we note that v = (3 — [n/ — l])/2. Hence, we see that increasing Trh, kmax/h, and Hi 
while decreasing Ng and is what we want to maximize nj. Clearly, Ng cannot be decreased 
beyond the minimal number of efolds Nmin that is necessary for a successful inflationary scenario 
(constraint 3) for a fixed Trh- This will be one of the strongest constraints for bounding nj — 1. 
Increasing Hi while keeping Ng (and Trh) fixed requires decreasing tg, since Ng Hiitg ^ki)- 
However, because Nmm also changes if Hi and tg changes, it is not possible to keep Ng fixed right 
at the constraint boundary. As Hi keeps increasing, it eventually runs into the tensor perturbation 
constraint 7. Also relevant for the case of low reheating temperatures is the fact that for sufficiently 
large Hi/TRH, we run into constraint 4. For each Trh, nj—l can be maximized through the Hi and 
tg variations subject to the constraints just described. 
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As Trh is increased towards the highest temperatures consistent with energy conservation, 
constraints 5 and 6 become relevant. Even though constraint 6 is a tiny bit stronger of a constraint, 
it is very similar in numerical value to constraint 5. This is fortunate because as described before, 
constraint 6 is imposed for computational convenience and constraint 5 arises from fundamental 
principle of energy conservation. In the parametric region where constraints 5 and 6 

compete, the reheating scenario is somewhat unrealistic in that the reheating time scale is very fast, 
taking the system away from the coherent oscillation perturbative reheating regime. However, to 
put a conservative upper bound, we account for this extreme parametric region as well. It is in this 
sense that the bound that we will obtain for the maximum n/ is reheating scenario independent. 

We next note that for lower reheating temperatures satisfying constraint 7 

rh bound^ 


(with Hi maximized to maximize nj—l). Hi has to be brought down when Tj^h is brought down to 
satisfy constraint 4: 

1 

2(«^-4)+l 

(40) 


Since we have saturated constraint 9, we see that 



% = 


2 



f 2 ^ V(v)| ^ ( Hi/{27t) \ 

y j \ xoiki) j 


^max 

ki 


\-v 


(41) 


increases as Hi is lowered. Hence, depending in particular on the numerical values of E/tmax 
kmax/ki, the required CO^ can exceed unity, violating constraint 8. 

Finally, if we choose kj/kmax ^ 10^^ and ^ 3 x 10^^, we can always satisfy constraint 

10 if n/ — 1 ~ 1. The current scale invariant isocurvature perturbation bound is given by 


A2(100fc,-) 

A2(100k,) 


< 3 X 10^2_ 


(42) 


If this scale invariant spectrum bound is assumed to bound the blue spectrum as well, we have 


Ek . ~ ^3 X 10-2 X 10-2 _ 10-2 


(43) 


for «/ — 1 ~ 1. Hence, we see that if we choose E^^^ < 1 and ki/kmax ^10 ^ , we can satisfy 
constraint 10. Therefore, we will focus on this parametric regime and ignore constraint 10. 
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Let us now find an explicit estimate of the largest value of n/. First, we consider the case of 

Trh > Tfiff (44) 


7 

5 / o \ ~ir~A 7 I 

r„^ = 22-^(A|j Mpr-^+^, 


2 X 10i°r3/4y:^GeV 

5^g*iTRH) ] 6 g X 10^5/4y:^GeV = 6, 

(45) 


for which constraint 7 becomes relevant. Saturating constraints 3 and 7 in the current approxima¬ 
tion scheme, we find 

Trh 


fVe~52.8 + xln- 


2 

— - In 


(46) 


3 IQiOGeV 3 IQiOGeV 
Although Trh appears here suggesting Trh should be minimized to maximize the n/ — 1 bound, 
the Trh dependence shown explicitly in Eq. (|3^ dominates. As constraints 5 and 6 are similar in 
magnitude, we use constraint 5 to maximize Trh for simplicity for this simplified analytic estimate. 
In other words, here we estimate 

10\ 1/4 


max Trh 


-A Mp 
g*J 

6.5 X 101^1/4 GeV , 


-Aj 
7t \l 2 ^ 


(47) 

(48) 


in which we have taken g* = 200 and A^(k,) = 2.4 x 10 From Eq. (|38|), we then find that 

55 


L7fi//=max7fiH 


< 

rs_/ 




^ ^ 55 


iXo)t, . R 


iv±p 




L 2 X 103 .3/4 

I jWor 710- N'/" 






jy 


1.2 X (l + ^ln 

■ .3/4 ■ 
. ^^max . 

) 


1-5.5X 10-3lnr-l.l x 10-2ln 

kmax —5 



(49) 


(50) 


in which we have taken g* and A|(k,) to have the same numerical values as above. We have also 
taken Xoi^k,) = Mp to satisfy constraint 12, we have approximated n/ — 1 1 on the right hand 

side, and we used 


1 2 / inr' 

Af' = 52.8 + -ln-J-In-- ^51.2( 1-. (51) 

^ 3 lOiOGeV 3 lOi^GeV V 205 ' 
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We see that there is only a very mild dependenee on the phenomenologieal parameterizations 
used above. Henee, if a blue isoeurvature speetral index nj >3 is measured, this eertainly eannot 
arise from a linear speetator with a time-independent mass. We will sharpen this estimate with a 
numerieal analysis in See.|^ 

We next eonsider the ease of 


/ 10 \ 


1/4 


M, 

but still with Hi saturating eonstraint 7: 


^6.5 X lO^V^/^GeV 


(52) 




1.06 


l-8x 10-3 lnr+8x 10-3 -1.2x 10-2ln[^ 10-5 

lOl'GeV L V 


X 


1 + 




+ Ain 


Trh 


44 lO'iGeV 


(53) 


where we have used 


isjest ^ 


AH A 1^1 

47.5-Inr-f - In—n-■ 

3 3 IQiiGeV 


(54) 


We ean eontinue to lower the temperature towards unless eonstraint 8 is saturated. Constraint 
8 is saturated before reaehing if «/ — 1 is smaller than the solution n‘j{T > Tn^) — 1 to 






("/-!) 


Tt 


3/2 


K^/rT 


3-(n/-i: 

2 


= 1 , 


(55) 


which for = 1, Wx/Vin = 10^Zo(^/t,■) = Mp.Vb = 10 is approximately 0.7. 

Let us now consider the lower reheating temperature Trh < Tn^ (and n/ — 1 > nf — 1), which 
means that we should set 


Hi ~ bound 


(56) 


in Eq. ([38]) instead of using constraint 7. We 





' 1.06(l-0.0231n[£,^„]+0.0381n^^^-7.6xl0-3ln[S87r]) ^ 

l+5.4xl0-3ln[58;r] 2.7xl0-3ln [f« 1.2x 10-2ln|^'='”“ lO-sJ 

1.2(l-0.0211n[£,„„]+0.0301n^^^-4.2xl0-3ln[58;r]) 

1+3.5X 10-3 ln[58;rl-rl.7x 10-3 In 1.3x lO-’ln 

kmax 2Q—5 

— \J 


At lower Trh, we have > 1 at 

( 4.6x10’ GeV ^/^ 

^max 


TdM2 = < 


, 105 J 

28 GeV eI'^ 

_ «max _ 

/max//njjn ) 5/^^ 

105 ) 


!-(«/- 2 ) 


8.5 + |ln[ ^™;^ff™ 


l-(n/-2) 


14+|ln 


^max / ^ min 


105 


riff = 5 
riff = 6. 


(57) 


(58) 
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We have taken the minimum Trh in this paper to be at 100 GeV to simplify the presentation. This 
means that (0^ < I constraint is more relevant for — 5 case than the — 6 case. 

We should also estimate the effect of constraint 11 on kmax^ 


which becomes 


kmax ^ ^ 

kmax 


h 


<10 


8 


(59) 


(60) 


with Xoihi) = Mp, nj—\ = 1.3, and Hi set at = 0.1. 

Finally, constraint 13 can be shown to be generically satisfied in the n/ — 1 and kmax/kj^jn 
region of interest. This will be discussed more in the numerical section below. 

3.2. Perturbative in Slow-roll Evolution 


In this subsection, we examine the effect of turning on £j^.. We will see its most important 
feature is to have nj maximized for \Hi — Hg\/Hg -C 1, making the Bessel spectral formula accurate. 

Instead of completely neglecting £k. during inflation in computing {{Xo) )t=tRH^ we can use 
linear perturbation theory in £k. to solve Eq. Q (for more details, see appendix [A|): 


{{Xo?)t>t^,tg ~ yZo (^fe,)exp [-(3 - 2Vk^)H,{tg - 1^.)] 


H{t) 

H{tg)_ 


= '* 16 v 3 - 3 ( 3 - 2 v ,)£, ) 

Fi = HiVkj{tg - tk.) - 1 

^ 32v4 + 6£k,Vk, {3 - 2Vk, [1 + H,{tg - 1 ^,) (3 + v^, [2 + (3 - 2Vk,)H,{tg - .)])]} + F 3 

F 2 = - - - - - 

^9-4v2 (i6v 3 -3e^,(3 -2v,,)(l +2Fi//,Vfc,(t, -t,,/ 

F3 = 9(3 - 2Vki)£k.Hi(tg - +2FiHiVki(te-^ki)], 

in which tm is the time after x field starts to oscillate and tg is the end of inflation: 

3 , , 
m = -H(t,n). 


This allows us to rewrite the analog of Eq. (381 as 




55 

mm 


^ ^ 55 


10 


(0^ min 


-( ) 

>/-l) Vl0 i°GeVy 


^55^^ 


iXo)l . rl(ni-l)R 


Mj 


(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 


( 67 ) 
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in which 




COxmmini-l) =max 


1 '\/l 


£h 






^min ) 5 




r,„(n/- 1) = 


m 






( 68 ) 

(69) 

(70) 


Nmin is given by Eq. (12) is given by Eq. (62), \//% is given by Eq. (jh I, and 

in 1 /~. z 

— -7(n/-l)(6-[n/-l]). 


(71) 


Comparing with Eq. (381, we see a eomplieated funetion r„j(n/ — 1) that depends on n/ — 1. Most 
of this eomplieated funetion aeeounts for the dependenee of the time evolution of 2^o(0- 

When aeeounting for e^. and eonstraint 2, we note that a given pair He and Ng can originate 
from two different values //, : 




H,= 


2 , H}Ne 


% 


M2a2 




(72) 


With HgNg Tt^MpA^, the hybrid inflation ease eorresponds to the minus sign braneh while the 
quadratie inflation ease eorresponds to the positive sign braneh. One ean also easily show that 
for the parametrie regime of interest, Hg never beeomes elose to zero even though there may be a 
worry from the form of Eq. ( [T0| ) that we may be unreasonably extrapolating the linear expansion 
of the slow-roll that is valid near On the other hand, the parametrie regions where the hybrid 
inflation and quadratie inflation branehes merge are sensitive to the branehpoint singularity there. 

The most important feature of turning on e^. is that sinee now (with eonstraint 2 imposed) 


An^A^pM'i /ff2 


(73) 


the minimization of Ng that is important for the extremization of nj — I (e.g. see Eq. ( [38| ) whieh in 
turn is related to eonstraint 9) gives a numerieal pressure in the non-linear extremization problem 
to make Hi elose to Hg. This favors a smaller (in turn favoring small //,) whieh eompetes with 
the pressure to extremize Hg (favoring a large //,) that arises from eonstraint 3. Henee, depending 


on the size of r, may not quite saturate eonstraint 7 as was done in the derivation of Eq. (50). 
This means that with Ej^. turned on, the sensitivity to the tensor-to-sealar ratio r entering eonstraint 
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7 is reduced for values of r that are “large”. As we will see in Sec. Q this makes the approximate 
spectral index Eq. ([^ more accurate for = 0.1. 

A figure of validity for the £]^. perturbations can be written as 

Cpert = £kN^, (74) 

since the background field evolution equation during inflation 

d^Zo + ^Hi{l-ek.Hi{t-tk.))dtZo + m^Zo = 0 (75) 

has a secular term £kHi{t — t^.), and this term is integrated over a time period of tg — tk^ Ng. 
Since Ng ~ 50, high Hi models where £, approaches the tensor-to-scalar ratio r bound have Cpgrt 
approaching unity, and hence they cannot be addressed reliably using this perturbative approach. 
In the next section, we will turn to a numerical analysis of this extremization problem, which will 
allow us to get a handle on situations such as these when perturbative methods fail. 


4. NUMERICAL RESULTS 


In this section, we perform a numerical analysis to find the largest«/ consistent with constraints 
1 through 13. The results of this analysis will show that even with an extremely optimistic experi¬ 
mental sensitivity of lO^^A^ on length scales as small as 10 kpc scales, the theoretical prediction 
from a constant mass isocurvature field scenario is that experiments will not measure spectral 
indices nj greater than 2.4. 

We begin with Fig. which shows the case in which ~ l^^max/^min = 10^; = 

10~\ 10^^; — 5,6; S = (8;r)^\ 10^^(8;r)^^}. The results show that the maximum temper¬ 


ature estimated in Eq. (48) agrees with the right end of each plot to better than 30% and the 


maximum n/ — 1 agrees with Eq. (50) to better than 5%. For the = 0.1 plot (the left plot), the 
reason why there is a drop of nj — I near ~ 5 x 10^^ GeV is due to constraint 5 (reheating 
energy conservation at time tg) pushing up //, as Trh is raised.^ This upward push of Hi is al¬ 


lowed because from the discussion around Eq. (73), constraint 7 may not be saturated depending 
on the size of r. This non-saturation is indeed the case for most of the r^, = O.l curve (which we 
have also checked directly numerically) and makes the approximate spectral index Eq. (1^ more 


This increases which in turn increases the split between //, and Hg. This then increases Ng, as can be seen in 


Eq. (73 I, under the assumption that the increase in the split is the most important effect. 
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Figure 1: The maximum measurable spectral index nj — I assuming = 1, = 10^; = 

10^\ 10^^; riff = 5,6; S = (87r)^\ 10^^(871)^^} is plotted as a function of Trh- (left) The bound on the 
tensor to scalar ratio ri, has been taken to be 10^'. The maximum spectral index with the amplitude pa¬ 
rameterized by {Fjtmax = 1) ^max/kmin = 10^} is around u/ = 2.11. For T < 10^ GeV, the lower solid curve 
corresponds to the reheating non-nonrenormalizable operator dimension of = 5 while the upper solid 
curve corresponds to the case of the non-nonrenormalizable decay operator dimension of = 6. For 
T > 10^ GeV, the two curves merge. The dashed curve corresponds to weakening the coefficient of the 
non-nonrenormalizable operator by a factor of 10. The vertical curve on the left portion of the boundary 
curves occur because the expansion rate there is too small in that parametric regime to produce measurable 
isocurvature perturbations (i.e. constraint 8). The dotted curve corresponds to evaluation of fiducial spec- 
fral index af kmax instead of The correction is small (except at the highest Trh where the dip occurs) 
because constraint 7 is not saturated for = 0.1. (right) Similar to the left plot except with a probably 
possible future bound of = 10^^. 

accurate. We see how the dotted eurve matehes the solid eurve exeept at the highest tempera¬ 
ture where the dip oeeurs, as we will discuss more below.^ For the = 10“^ ease, constraint 
7 does saturate at the highest allowed reheating temperature, whieh means that no upward push 
of Hi ever arises from constraint 5 for these highest temperatures. The maximum nj — \ for this 
~ ^5 ^max/^min = 10^} experimental seenario is about 1.1. Any measurements of CDM- 
photon blue isoeurvature with a speetral index larger than «/ = 2.1 with an amplitude larger than 

^ The bottom of the dip is where the mismatch of the accurate dotted curve and the approximate solid curve is the 
largest. This does not affect our main result since it does not correspond to globally the largest spectral index. 
Furthermore, this reheating sliver is where the reheating scenario is least realistic and has been considered only to 
give a conservative bound on «/. 
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equal to = 1, ^max/^min = 10^} imply the responsible dynamical degree of freedom during 

inflation cannot be a constant mass linear spectator field. 

Let us now consider some of the other features of these results. For all but one of the curves 


shown in Fig. when Trh from above to below Tng while nj —I > — 1 (see Eq. (551 for the 


definition), there is a break in the bound curve as expected from Eq. (56) encoding the minimal 
reheating constraint 4. The break in the curve does not exist for the case of = 10“^, — 6} 

because in that case «/ — 1 reaches rfj{T > Tn^=e) — 1, in which 


7;^=6(^fc= 10~^5= (8;r)”i) ^ 120GeV, 


(76) 


which means that the dark matter constraint 8 is saturated without saturating the reheating con¬ 
straint 4. All of the curves terminate at a certain lower endpoint of reheating temperature because 
of the dark matter constraint 8, which simply states that the expansion rate in that parameter regime 
is too small to produce measurable isocurvature perturbations. We note that there is no vertical 
line plotted for the right hand side of the curves in Eig. [^(unlike the left vertical line) because we 
did not want to obscure the drop in n/ — 1 for high Trh for ri, = 0.1. 

Eig. shows the case with kmax/^min = 10^; = 10“^ 10“^; = 5,6; 5 = 

(8;r)^\ 10^^(8;r)^^}. Decreasing ^ means increasing the experimental sensitivity 

(i.e., A^/A| is resolved to 10^^ instead of order unity - an extremely optimistic view of the 
foreseeable future that is chosen to illustrate the insensitivity of the bound to experimental pre¬ 
cision). This changes the measurable maximal spectral index logarithmically to about nj = 2.25 
(from nj = 2.1 when = 1). Hence, although increasing experimental sensitivity changes the 
measurable blue spectral index, the logarithmic nature of the increase makes these numbers exper¬ 
imentally meaningful for at least a many decades time scale. As before, the maximum temperature 


estimated in Eq. (481 agrees with the right end of the plot to better than 30% and the maximum 


ni — I agrees with Eq. (50) to better than 10%. Eor the = 0.1 plot (left plot), the reason why 
there is a drop of nj — I near Trh ~ 5 x 10^^ GeV is the same reason as in the explanation for 
Eig.[^ Note that unlike in Eig.[^ the bounds for nff = 6 end at Trh = 10^ GeV because we simply 
truncated the plot there (and not because > 1 there). 

Einally, to be extremely optimistic regarding short distance scale probes of cosmology, in Eig.|^ 
we consider the nj—l bound with an experimental probe length scale of kmax/kmin = 10^ (i.e. kmax 
is at the scale of 10 kpc) with the other parameters set at {Eyt^ax ~ 10~^, = 10~\ 10~^; = 

5,6; 5= (8;r)^\ 10^^(8;r)^^}. The maximum spectral index increases as expected in a mild 
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Figure 2: The maximum measurable spectral index nj — \ as a function of Trh G [100,5 x 10^5] GeV 
assuming an experimental sensitivity of = 10^^ corresponding to resolving A^/A^ to 0(10^^%) at 
about 1 Mpc length scale. The rest of the parameters are set at {/cmax/^min = 10^; = 10^', 10^^; riff = 

5,6; 5 = (87r)^\ 10^^(87r)^^}. (left) The bound on the tensor-to-scalar ratio r has been taken to be 10^^ 
The lower solid curve for Trh ^10^ GeV corresponds to the reheating non-nonrenormalizable operator 
dimension of = 5 just as in Fig.[^ The dotted curve corresponds to evaluation of fiducial spectral index 
at kmax instead of kmin- As in Fig. the correction is small because constraint 7 is not saturated even 
for = 0.1. As expected from increasing the experimental resolution by 10^, the maximum measurable 
spectral index has only gone up mildly to nj = 2.25 (from nj = 2.12). Note that unlike in Fig.[T] the bounds 
for Jiff = 6 end at Trh = 10^ GeV because we simply truncated the plot there (and not because ( 0 ^> \ 
there). The dashed curve corresponds to weakening the coefficient of the non-nonrenormalizable operator 
by a factor of 10 just as in Fig.[^ (right) Similar to the left plot except that rj, has been set to 10^^. 


manner to r/ = 2.35 (from «/ = 2.25 with kmax/^min = 10^)- Note that this ^max/^min lies near the 


edge of constraint 11 in accordance with Eq. (60). 

Also, constraint 13 can be shown to be generically satisfied in the «/ — 1 and ^max/^niin region 
of our interest. For example, Fig. |^shows as a function of r/ — 1 defined according to 

Eq. ([^ for the parametric choices of = 0.1 and ^max/^min = 10^ (which is the most constrained 
among the scenarios we are interested in). We see that since we have considered only r/ — 1 < 1.6, 
constraint 13 will be satisfied. 

Note that for the numerical computations discussed thus far, only the background fields are 
evolved fully numerically to determine (0^ while analytic approximations relying on H{t) being 
constant have been used to compute A^/( 0 ^ in accordance with 117311 . For r ^ 0.1 we have 1, 
and H{t) evolution does not present much of a correction. However, for = 0.1 in the plots 
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Figure 3: Similar to the left Fig. [^except with A:max/^min = 10^ (i.e. /cmax is at the scale of 10 kpc). The 
bound on the maximum spectral index «/ is only logarithmically sensitive to kmax/^min as it is now about 
2.35 instead of 2.25. The features of the plots are explained as in previous figures. As before, the dotted 
curve corresponds to the evaluation of the spectral index with the fiducial k value of kmax instead of k mi n. 



Figure 4: is plotted as a funcfion of u/ — 1 defined according fo Eq. ^ for fhe severesf paramefric 

choices of = 0.1 and kmax/^min = 10^. This shows consfrainf 13 is salisfied for u/ — 1 < 1.6. 

above, there may be a worry that the numerieal eomputation of would deviate signifieantly 

from the approximations. One symptom of the analytie mode funetions destroying the aeeuraey 
of ean be tested by eomparing the answers for two different fidueial values of ko. 

The parametrie speetral indiees n/ — 1 shown in Figs. through [^(exeept for the dotted eurves) 
eorrespond to the ^ approximate parameterization with ko ehosen at ko = ki whieh is 

the longest observable wave veetor. For Hi eorresponding to saturating eonstraint 7 with r < 
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10 this is a good parameterization: i.e., in Figs. through plots with = 10 ^ ean be 
taken to be aeeurate to better than 1%. However, if Hi saturates the limit of eonstraint 7 with 
r 10^ ^ H{t) would evolve nontrivially during inflation. In that ease, there is a worry as to 
whether the {k/kiY‘^^ parameterization is inaeeurate for the = 0.1 eases. For example, if we 
saturate eonstraint 7 with = 10“^ a more aeeurate approximation of the observed speetrum near 
k max should have the fidueial value Uq = k^ax (at the expense of eomputing Xoiho) numerieally). 
Fortunately, we find numerieally that eonstraint 7 is never saturated even with r/, = 0.1 beeause 
of the effeets diseussed in Eq. ( |7^ . The aeeuraey of the analytie speetrum ealeulation ean also 
be seen in the dotted eurves of Figs, [^through[^whieh were eomputed numerieally.^®By explieit 
eomputation, we have eheeked that nj — I eomputed with mode funetion evolution evolved fully 
numerieally matehes the r/ — 1 eomputed through the Bessel funetion with ko shifted to kmax (and 
eomputed numerieally) to better than a few pereent. 

Henee, we eonelude from Fig.|^that any measurement of r/ > 2.4 for CDM-photon isoeurva- 
ture perturbations in the foreseeable future indieates the responsible dynamieal degree of freedom 
during inflation eannot be a eonstant mass linear speetator field. 


5. MODELS: WHAT HAPPENS WITH A DYNAMICAL MASS? 

In [l73l . it was shown that speetral indiees as large as r/ = 3.8 (but not r/ = 4) ean be aehieved 
in the eontext of a dynamieal VEV breaking the Peeeei-Quinn (PQ) symmetry. This is of interest 
beeause r/ = 3 is eonsidered to be observable for example by the Square Kilometer Array iTT^ . 
Here, we diseuss why a time-dependent mass during inflation ean evade the bound diseussed 
around Fig. Suppose the mass of the field x responsible for the linear speetator isoeurvature 
makes a transition at time tc from value m to zero. Aeeording to eorollary 2 of iTT^ . the modes 
k < kc that leave the horizon earlier than the time of the mass transition still have the form of 
Eq. Q, where if the slow-roll evolution is negleeted, the eritieal wave veetor is given by 

~ ^min6Xp[Ac], (27) 

in whieh 

(78) 

The agreement between the dotted curve and the solid curve exists except at the highest Trh dipping sliver which 
does not correspond to the globally maximum «/, as we discussed in footnotel9| 
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is the number of efolds from the beginning of inflation. For kc ~ the number of efolds 

for whieh this oeeurs must be at least A^c ~ 16 efolds. All of these modes are governed by massive 
sealar field quantum fluetuations giving a blue speetrum. In addition, the baekground field Xo{^) 
dilutes as 

nr — 1 

Xoit) exp[-^^c], (79) 

which dilutes the total isocurvature by an important factor: 

A2ocexp[-2(R/-l)Ac,] (80) 


which is analogous to Eq. ( [^ . The field theory up to this point behaves just as in the constant 
mass scenarios we have been discussing. 

However, after the mass transition to masslessness completes, the background field Xoi^) be¬ 
haves as a constant massless field until the end of inflation. Hence, compared to the constant mass 
case, the isocurvature perturbations receive a boost of 


A^(time dependent mass) 


exp[2 (n/- 1) {Ne-Nc)] 


(81) 


A^ (constant mass) 

in which Ne is the total number of efolds as usual. Since Ng ~ 50, the enhancement for Ac ~ 16 
scenario is enormous. This is the intuitive explanation with which time-dependent mass situations 
can evade the blue spectral index bounds for the time-independent mass situation that has been the 
main focus of this paper. One observational signature of the mass transition iT/^ is the existence 
of a flat isocurvature spectrum (for k > kg) in addition to the blue spectrum {k < kg). On the other 
hand, if there is a limited k-range accessible experimentally, it may not be easy to observe the 
break in the spectrum. 

A natural question is then what class of models naturally produce these time dependent masses. 


Note that the crucial ingredient in being able to generate the large enhancement Eq. (811 is the 
transition from m///,- ~ 0(1) to m///; ^ 1. If //,■ is the natural minimum energy scale for the 
masses of the scalar dynamical degrees of freedom (as is the case for example in supergravity 
models), then a symmetry needs to naturally lead to m///,- -C 1. Hence, one crucial ingredient for 
natural isocurvature models with n/ larger than the bound presented for the constant mass case is a 
symmetry X protecting the X mass from Hubble scale corrections to its mass. A second ingredient 
is a temporary (but lasting many efolds) mass generation mechanism. This second ingredient is 
necessary to generate the blue spectrum. 
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In the supersymmetric axion scenario of 11741 . the symmetry X is the Peccei-Quinn (PQ) sym¬ 
metry non-linearly realized as a shift symmetry of the axion field. The PQ symmetry breaking 
fields are displaced from the minimum of the effective potential during inflation (in a way in 
which PQ symmetry is always broken) such that the coset symmetry X is actually broken by 
through the kinetic structure of the axion: i.e., the Nambu-Goldstone theorem does not apply be¬ 
cause the system is not in vacuum. As fields roll toward the vacuum (where the PQ breaking 
persists), the axions behave as a massive field with mass of the order of Hi due to the supergravity 
structure of the Kahler potential. After reaches the vacuum and the kinetic energy dilutes to 
the point of ^ X is restored axions become massless, up to the small explicit PQ 

breaking contribution. 

Although it is possible to tune parameters and initial conditions to obtain almost flat potentials, 
the Nambu-Goldstone models with out-of-equilibrium symmetry-breaking time-dependent VEVs 
seem to be the simplest natural model. From this perspective, any experiment measuring a CDM- 
photon isocurvature perturbations with n/ > 2.4 may be finding evidence for a dynamical degree 
of freedom during inflation that has a coset shift symmetry. 


6. CONCLUSIONS 

We have considered a constant mass spectator linear isocurvature degree of freedom during 
inflation and answered the question of what is the largest measurable blue spectral index that can 
be produced via such a mechanism. We have shown that the largest measurable spectral index is 
less than 2.4 in the foreseeable future with only logarithmic sensitivity to experimental precision 
characterized by {F/tmaxi ^max/^min} and experimental constraints such as the tensor-to-scalar ratio 
r^. This means that any future measurements of the isocurvature spectral index above this bound 
would give weight to the hypothesis that there is a spectator field with a time-dependent mass 
during inflation. 

We have also considered how for reheating temperatures much smaller than the maximum 
allowed by tensor perturbation bound, the maximum observable spectral index decreases. This 
would be relevant if there were specific inflationary models under consideration with a fixed re¬ 
heating scenario or model-dependent phenomenological bounds on the reheating temperature such 
as those that arise from cosmologically dangerous gravitinos. For part of this smaller reheating 
temperature dependent bound, we have used the assumption that there is at least a gravitationally 
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suppressed non-nonrenormali / ahle operator of dimension 5 or 6 that ean eontribute to reheating. 
This assumption sets a bound on the maximum separation between the reheating temperature and 
the expansion rate at the end of inflation in eertain eases. 

One has to keep in mind that the maximum derived in this paper has some obvious eaveats. 
First, sinee we have only eonsidered linear speetator seenarios, we have not examined what the 
maximum blue speetral index would be if we allowed to be of order Xq . Sinee we have imposed 

Xq > H/(In) and dx is at most of order H/{In), one might think that the eurrent estimate will 
stand even after ineluding the 8x^ Xq seenarios. On the other hand, the «/ — 1 of quadratie 
isoeurvature seenarios (i.e. seenarios in whieh the isoeurvature perturbations are proportional to 
twiee that of the linear speetator seenario 11100111021 . However, a preliminary 
investigation shows that this faetor of 2 in the power only gives an enhaneement of the form 

max n/ - 1 oc-i (82) 

l-^ln[kn,axAniin] 

multiplying a diffieult to eompute suppression (originating from the quantum nature of the parti- 
ele produetion in eontrast with the elassieal VEV displaeements of the linear speetator seenario), 
resulting in a similar maximum speetral bound at best. However, given that the dependenee of the 
relie density and the speetral amplitude with n/ — 1 is somewhat eomplieated due to their depen¬ 
denee on the long time mode evolution lllOOl . it would be worthwhile eonfirming the quadratie 
isoeurvature estimate more earefully. 

Another eaveat is that we have assumed a “standard” slow-roll, effeetively single-field infla¬ 
tionary seenario with only one reheating period. Most non-minimal extensions will dilute the 
VEV energy density leading to a smaller upper bound. In that sense, most of the non-minimal 
extensions are not likely to ehange this general pieture. Even in the situation in whieh Xo makes a 
phase transition after inflation (e.g. Xq goes from vi to V 2 ) sueh that CO^ (now proportional to v^) 
is generated after inflation (thereby evading the inflationary dilution), sinee it is really dx that is 
diluting during inflation (even though we have been rewriting it as C0;(; 5;( / |;(o ~ ^ 1 1 being eonstant 
during inflation), this does not help us to evade the bound. 

Einally, we have assumed a sampling of inflationary spaee eharaeterized by {£i,Hi,te}, while 
there are infinitely more ways to tune the inflationary models. On the other hand, even the addition 
of £; (versus a non-evolving seenario of H{t) during inflation) produeed only about a 10% ehange 
in n/ — 1. Henee, we believe this limitation of sampling is not severely restrietive. 

It is indeed intriguing that future eosmologieal inhomogeneity measurements of r/ > 2.4 may 
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uncover the following new features of dark matter: (i) dark matter had to have a time dependence 
in its mass in its evolution history in the context of an inflationary universe, and (ii) dark matter 
mass was of order of the expansion rate during inflation. From our current model building tool-kit, 
arguably the most appealing picture that would emerge is that the dark matter is a field possessing 
a fundamental shift symmetry just like the axion. 
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Appendix A: Background Solution 


Eor this section, we set the time at which the observable longest wavelength mode leaves the 
horizon to be time = 0. We can model a very large class of slow-roll inflationary models 
with the Hubble expansion rate function parameterized (with three constants {eytpi7i,L}, where L 
approximately replaces rjy in the usual slow-roll parameterization scheme) as 


H 


Hjjl-ek.Hjt,) 
l + lit-teWl-Hfy.te) 


0<t <te 

te <t < tRH 


(Al) 


in which At = t — ti and tRH is the time of reheating. This ansatz accurately (at the order of 10% 
level) both quadratic inflation and hybrid inflation. Note also that as long as the number of efolds 
is fewer than 


N„ 


2et 


AK^MjA^iki) 

hF 


(A2) 


the quantity H will never go negative. 

After the end of inflation, the solution of the field evolution equation 


takes the simple form 


Xo{t) + 3HXo{t)+mlxo{t) = 0 


Xoit) 


eiCOs{m),At) + e2sin{mAt) 

l-lHiAt{ek,Hite-l) 


(A3) 


(A4) 
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where At = t — tg- 

We eould in prineiple solve the equation of motion (Eq. ( |A3[ )) exaetly in this elass of models in 
terms of hypergeometrie funetions and Hermite polynomials 

—1 


Xo = CiH „,2 ( 


2£k, 


+ \ I +C2 iFi 


ietflr r \V let, Vz*'*'"'' 


(A5) 


However, beeause is small, these speeial funetions must be evaluated in exponentially large and 
small numerieal regions and added together. Sueh a route seems numerieally unstable, in addition 
to being opaque. In praetiee, it is easier to handle numerieally the solution to the equation of 
motion subjeet to the boundary eondition 

Zo (0) = - Q - (0), (A6) 

whieh embodies the assumptions that the speetral index is of order unity and the field is rolling in 
a slow-roll fashion, initially. 

We ean mateh the solution before and after the end of inflation to write the solution after the 
end of inflation as 


Xo{t)=Ki 


m 


Wte) 

Ki = A/^;i:o(0)exp 

in whieh K 2 is a phase. The amplitude is given by 


cos{mAt + K2) 

-ljp>-2Vki)Hite 


(A7) 

(A8) 


Zo(0) \ 


1 ^kiHite~\~ 3Zo(^e)/(^iZo(^e))] 

^ I ^ A ^ 


1 - H 


(A9) 


in whieh we note that Xoih) exp [^(3 — 2v,)//;te] is the initial value ;to(0) for ^ki = 0- Henee, it 
is more eonvenient numerieally to solve for Xoi^) exp [^(3 — 2v,)//;t] than Xoi^)- The exponential 
suppression of Xq (4) exp ~ 2Vi)Hite\ / Xq (0) still oeeurs when 9 — An? jU ^/(1 ~ £k Hi tef<0. 
In this notation, the dark matter fraetion co^ is 

Kj m2 Trh 


% — 




-R. 


(AlO) 


eq 


where R is defined in Eq. 25 
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